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Abstract 

Expressions  for  the  radiation  patterns  of  a  slit  in  an  infinite  ground 
plane  covered  by  a  uniform  layer  of  plasma  are  obtained.  The  slit  is  con¬ 
sidered  to  be  infinitely  long.  Two  polarizations  are  considered.  In  one 
case,  the  electric  field  is  polarized  across  the  gap,  and  in  the  other,  the 
electric  field  is  polarized  along  the  slit. 

The  effects  of  anisotropic  plasmas  on  the  radiation  pattern  are  also 
studied.  However,  only  the  non-gyrotropic  plasmas  are  considered. 

The  solutions  to  Maxwell's  equations  are  expressed  in  terms  of  spec¬ 
trum  of  plane  waves  in  the  rectangular  coordinate  system.  Fourier  trans¬ 
form  relations  are  used  to  match  the  boundary  conditions.  The*solutions 
are  obtained  in  the  integral  form,  and  the  saddle-point  integration  method 
is  used  to  obtain  the  far-field  asymptotic  expression  of  the  solutions. 

The  normalized  power  patterns  are  plotted  for  the  various  cases 
studied,  and  the  results  discussed. 
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Radiation  Pattern  of  a  Slit  in  a  Ground 
Plane  Covered  by  c  Plasma  Layer 


1.  INTRODUCTION 

When  a  supersonic  missile  re-enters  the  earth's  atmosphere,  it  becomes 
surrounded  by  a  layer  of  ionized  gas.  The  presence  of  the  ionized  gas  around 
the  re-entry  vehicle  alters  the  performance  of  the  vehicle's  radiative  system. 
The  effects  of  the  plasma  layer  on  the  attenuation  and  the  shape  of  the  antenna 
pattern  are  difficult  to  predict  because  the  plasma  layer  is  highly  nonuniform, 
and  the  geometry  of  the  radiating  structure  is  not  always  simple. 

The  object  of  this  report  is  to  study  the  attenuation  and  the  changes  in 
the  far-field  radiation  pattern  due  to  a  plasma  layer  for  a  simple  radiative 
system,  namely,  a  slit  in  an  infinite  ground  plane  covered  by  a  uniform  layer 
of  plasma.  The  geometry  is  a  good  approximation  to  either  a  horn  or  a  slot 
mounted  flush  with  the  body  of  the  re-entry  vehicle. 

Assuming  that  the  slit  is  infinitely  long,  the  problem  is  reduced  to  a  two- 
dimensional  one.  The  plasma  layer  is  assumed  to  be  uniform  and  linear. 

Two  different  polarizations  are  considered.  In  Section  2,  the  case  of  the 

electric  field  polarized  across  the  gap  is  investigated.  A  limiting  case  of  this 

2 

particular  polarization  has  been  investigated  by  Newstein  and  Lurye.  The 
power  patterns  are  plotted  and  the  features  of  the  patterns  are  explained  in 
terms  of  simple  models. 

In  Section  3,  the  case  of  the  electric  field  polarized  parallel  to  the  slit  is 
studied.  Radiation  patterns  with  anisotropic  but  non-gyrotropic  plasma  layer 
are  investigated  in  Section  4. 


(Received  for  publication,  31  October  1962) 
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1.1  Method  of  Solution 

To  solve  the  problem,  the  solution  of  Maxwell's  equations  is  expressed 
as  an  infinite  space  spectrum  of  plane  waves  in  rectangular  coordinates. 
Fourier  transform  relations  are  used  to  match  the  boundary  condition  at  the 
ground  plane.  By  applying  boundary  conditions  at  the  other  interface,  the 
solution  of  the  field  is  expressed  as  a  Fourier  integral.  Finally,  to  obtain 
the  far-field  radiation  pattern,  the  saddle-point  integration  method  is  used 
to  evaluate  the  integral. 

3,  8,  9 

1.  2  Effective  Permittivity  of  a  Plasma  Medium 

In  this  section,  the  effective  dielectric  permittivity  for  the  simplest 
model  of  plasma  medium  is  derived,  and  procedures  to  obtain  the  permit¬ 
tivity  of  more  complicated  models  are  indicated. 

Plasma  is  an  ionized  gas  which  is  electrically  neutral.  Since  the  ions 
are  much  heavier  than  the  electrons,  the  velocity  of  ions  due  to  the  applied 
fields  is  negligible  when  compared  with  the  velocity  of  the  electrons  in  this 
model.  The  effect  of  collision  is  neglected  at  first. 

Assuming  a  time  dependence  of  the  form  e^'^*'.  Maxwell's  equations  in 
complex  form  are  given  by: 


V  xE  =  -pp^  n 

(1.  la) 

V  X  H  =  E  +  j 

(1.  lb) 

J  =  pv 

(1.  Ic) 

where  p  is  the  electron  charge  density  and  v  is  the  velocity  of  the  electrons. 

It  shall  be  assumed  that  the  medium  being  dealt  with  is  homogeneous  and  that 
the  applied  electromagnetic  fields  are  not  so  large  as  to  change  appreciably 
the  homogeneity  of  the  medium.  In  other  words, 
p  ~  Pq  =  constant. 

The  equation  of  motion  for  the  electrons  is  given  by 

m^  =  e  (e  +  V  Xpq  h)  (1.  2) 

where  e  =  election  charge  and  m  =  electron  mass.  With  the  initial  assumption 
that  there  is  no  magnetostatic  field  present,  the  quantity  v  Xp^H  becomes 
second  order  and  is  neglected.  The  equation  of  m.ution  in  complex  form  reduces 
to 

jwmv  =  eE 

-  e  E 
jwm 


V 


(1.  3) 
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Substituting  this  relation  into  Eqs.  (1.  lb)  and  (1.  Ic), 

E 


V  X  H  =  +  — — 


0  jajm 
^0® 


=  ]ae 


0 


1  - 


2 

u.  “me 


0  J 


But  Pq  =  N^e  where  is  the  electron  density.  Therefore, 


V  X  H  =  ia  e 


0 


1  - 


"0 


This  gives  the  effective  dielectric  constant  of 


P  0 


,  ,  N  e^ 

where  2.  0 

u)  =  - 

p  me 


^0^ 


2 


a>  ^ 


0 


(1.  4) 


(1.  5) 


(1.  6) 


Since  the  velocity  of  the  electrons  is  90°  out  of  phase  with  the  applied 
electric  field,  there  is  no  loss  (  see  Eq.  (1.  3)]  .  When  the  effect  of  collisions 
is  appreciable,  the  velocity  of  the  electrons  is  no  longer  90°  out  of  phase  with 
the  applied  field.  The  effect  of  collisions  is  included  in  the  equation  of  motion 
as  a  damping  term  given  below: 


d7 

dt 


+  m  rv  =  e  E 


(1.  7) 


where  v  is  the  collision  frequency.  For  such  a  system,  the  permittivity  is 
^ven  by 


2  -I 


2^  2 
V  +(j.' 


,  ■  ^  I 

+  3-2 


(1.8) 


Macroscopically,  the  important  difference  between  an  ordinary  dielectric 
and  a  plasma  medium  is  that  the  real  part  of  the  permittivity  for  the  plasma 
medium  is  less  than  e^.  Actually,  the  real  part  of  e  can  talte  on  negative 
values. 

If  a  magnetostatic  field  is  introduced,  then  the  equation  of  motion  for  the 
electrons  becomes 


m^  +  mrv  =  e  j^E  +  v  X  Bq  j  .  (1.  9) 

Such  an  equation  of  motion  gives  rise  to  a  tensor  permittivity  of  the  following 
form  (in  the  rectangular  coordinate  system): 
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The  magnetostatic  field  is  assumed  to  be  in  the  z  direction  in  a  right-hand 
coordinate  system.  The  quantity  is  the  gyrofrequency  and  is  defined  as: 

b  m 


1.  3  Plane  Waves 

In  the  analysis  of  the  radiation  problem,  the  fields  are  expanded  in  terms 

of  plane  waves.  A  brief  review  of  plane  waves  in  semi -infinite  medium  is 

i  wt 

given  here.  In  the  subsequent  analysis,  the  time  dependence  e  is  assumed. 

From  Maxwell's  equations,  the  following  equation,  known  as  Helmholtz 
equation,  can  be  obtained  for  the  electric  field: 

[  E  =  0  , 

where  for  the  rectangular  coordinate  system 


In  general,  the  solution  to  such  an  equation  in  rectangular  coordinate  system 
is  given  by: 


E(x,  y,  z)  ■  \  z  , 


where 


E,  ^  “  component  of  electric  field 

IX,  y j  Ti) 

Aq,  r^,  r.,,  and  Tg  are  complex  constants 

,  2  ,  ,  2  ,2  2 
k  +k  +k  =  ill  ue 
X  y  z 


(1.  11) 


Consider  the  geometry  shown  in  Figure  1.  Assume  that  a  uniform  plane  wave 
impinges  upon  the  semi -infinite  isotropic  dielectric  from  a  direction  making  an 
angle  0  with  z  axis,  and  that  the  electric  field  is  polarized  in  the  y  direction. 

With  these  simplifications,  the  solution  in  the  free  space  region  reduces  to 


(1.  12) 


Figure  1. 


The  first  term  in  Eq. 
wall  and  the  second  te 
eauations  the  maeneti 


6 


and 


BE 


H  — 
z 


= 


9x 


-k 

_ X 

OJfJp 


A’ 


-jk  X 


.  (1.  14) 


For  this  particular  polarization,  all  other  field  components  are  zero. 

Inside  the  dielectric  medium  there  is  no  reflected  wave.  Thus,  the 
expression  for  the  electric  field  is  given  by 


E 


T  e 


-j(k'^x+k'  ^z) 


where 

T  =  constant 


t  I 

k  =  k  sin  0 

X 


I  t 

k  =  k  cos  0 
z 


2 

oj  pe 


(1.  15) 


Again  the  magnetic  field  can  be  obtained  from  Maxwell's  equations.  From 
the  boundary  conditions,  the  value  of  F  and  T  can  be  obtained,  as  well  as  the 

I 

relation  between  0  and  0  . 

For  a  more  detailed  analysis  of  plane  waves  see  Reference  4. 

In  this  problem,  the  solution  to  Maxwell's  equations  will  be  expressed 
in  terms  of  plane  waves  like  those  of  Eqs.  (1.  12)  and  (1.  15).  ^ 

In  the  anisotropic  case,  there  is  a  different  relation  between  k  ^  and  k 
A  discussion  of  plane  waves  in  an  anisotropic  medium  is  given  in  Section  4. 
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2.  TRANSVERSE  MAGNETIC  WAVE 
2. 1  Solution  in  Integral  Form 

The  geometry  for  the  first  problem  considered  is  shown  in  Figure  2. 

At  the  gap,  the  electric  field  is  polarized  in  the  x  direction  and  the  amplitude 
distribution  along  the  gap  is  assumed  to  be  constant,  as  illustrated  in  Figure 

3.  The  magnetostatic  field  is  assumed  to  be  zero. 

A  very  similar  problem,  namely  that  of  a  magnetic  line  source  in  the 

ground  plane  covered  by  a  layer  of  plasma,  has  been  solved  by  Newstein 
2 

and  Lurye.  “  The  magnetic  line  source  in  a  ground  plane  is  a  limiting  case 

of  a  slit  in  a  ground  plane  with  the  slit  becoming  infinitesimally  narrow. 

2 

Contrary  to  Newstein  and  Lurye's  method  of  expanding  the  solution 
to  Maxwell's  equations  in  terms  of  cylindrical  waves,  the  solution  in  this 
report  is  expanded  in  terms  of  plane  waves  in  rectangular  coordinates. 
However,  the  results  are  the  same. 

It  is  clear  from  the  boundary  conditions  that  for  the  particular  geometry 
and  polarization  described  in  Figures  2  and  3  that 

E  =  H  =  H  =0 
y  X  z 

I  I  I 

E  =  H  =  H  =0 
y  X  z 

where  the  primed  quantities  represent  the  electric  and  magnetic  fields  within 
the  plasma  and  the  unprimed  quantities  represent  the  fields  in  free  space. 


FREE  SPACE 

*0  Mo 


Figure  2.  Radiative  Geometry  for  the  TM  Wave  Case 
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Figure  3.  Electric  Field  at  z  =  0  for  the  TM  Wave  Case 


In  the  free  space  region  z  >  1,  the  following  solution  is  assumed  for  the 
magnetic  field: 

noo  -jtkx+kz) 

H  (x,  z)  =  \  aCx  )  e  dk  , 

y  .Loc  X 

'"here  ^  -  k^2  .  (2.  1) 

a(k^)  =  amplitude  function  independent  of  x  and  z.  Thus,  the  solution 
is  expressed  as  a  spectrum  of  outward-going  plane  waves.  If  a  plane  wave 
satisfies  Maxwell's  equations,  then  certainly  the  sum,  or  more  precisely, 
the  integral  of  such  plane  waves  satisfies  Maxwell's  equations. 

It  should  be  noted  here  that  along  the  path  of  integration,  as  k  goes 
from  minus  infinity  to  plus  infinity,  the  quantity  k^  takes  on  both  real  and 
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imaginai'y  values.  The  sign  of  must  be  chosen  in  such  a  way  that  the 
imaginary  part  of  is  always  negative.  If  the  positive  values  are  taken 
for  the  imaginary  part,  the  integral  is  expre.ssed  in  terms  of  increasing 
exponentials,  and  the  integral  will  not  converge.  Physically  such  a  result 
is  not  expected. 


From  Maxwell's  equations  result  the  following  expressions  for  the 
electric  field: 


E  (x,  z) 

X  jtae. 


3H 
_ V 

az 


«00  , 


a(k  )  e 

X 


-j(k  X  +  k  z) 

X  z 


dk 


(2.  2) 


E  (x,  z) 
z 


_1 _ 


k^ 

9x 


-1  (k  X  +  k  z) 
\  X  z 

a(k  )  e 

X 


dk 


X 


(2.  3) 


Inside  the  plasma  medium,  the  solution  for  the  magnetic  field  is  repre¬ 
sented  by  a  spectrum  of  outward-going  waves  plus  reflected  waves  as  follows; 

t  !  t 

-ik  I 


1  “jk  Z 

1  jk  „z' 

i.  x)  e 

+  c(k  e 

dk 


(2.  4) 


where 


k  =  (k  “-k'  )‘ 

z  p  X 


k  = 

P  Op 


The  propagation  constant  k  is  in  general  a  complex  quantity.  The 
quantities  b(k  ^^1  and  c(k  j.)  give  the  complex  amplitudes  of  outgoing  and  re¬ 
flected  waves. 

The  corresponding  electric  field  is  obtained  from  Maxwell's  equations; 

-1  ’"'y 

P 


'  aje 
•  -“=  D  L 


-jk  z  I 

b(k  e  ^  -  c(k  e 


jk  ^z 


-jk  X  , 
j  ^  dk^.(2.5) 


,  jk' 
+  c(k  )  e 

X 


-jk  X  , 

e  ^  dk  (2.  6) 


T  I 

There  are  three  unknowns  in  this  problem:  a(k^),  b(k  ^),  and  c{k  ^). 
Tb  evaluate  them,  it  is  necessary  to  have  three  equations  relating  these 
quantities.  The  boundary  conditions  will  provide  the  necessary  relations. 
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Consider  the  boundary  condition  at  the  ground  plane.  Here  the  ampli- 

I 

tude  of  the  tangential  component  of  electric  field  h  as  function  of  x  is 
as  shown  in  Figure  3. 

I 

Equation  (2.  5)  for  E  ^  reduces  to  the  following  at  z  =  0; 


E  (x 

X 


poo  1, 

-  01=  \  77 


k/  —  00 CJ  £ 


b(k  J  -  c(k 


-jk  X 
e  dk 


(2.  7) 


Equation  (2.  7)  is  a  familial-  expression  for  the  Fourier  transform  between  the 

r 

k  ^  domain  and  the  x  domain.  It  can  be  seen  readily  that  the  quantity 


(j  e 


b(k  ^)  -  C(k 


is  the  Fourier  transform  of  the  electric  field  amplitude  distribution.  E  ^{x,  0). 
The  Fourier  transform  of  a  rectangular  pulse  is  a  simple  one,  and  thus,  the 
transform  will  be  equated  to  the  above  expression. 


k 

_ ^ 

we 

P 


)  -  c(k 

X 


2’rEQa  sin  k 

k’  ^ 
^  X  2 


(2.  8) 


Equation  (2.  8)  gives  one  of  the  three  required  relations,  and  the  other  two 
relations  are  obtained  at  the  interface  between  the  plasma  and  the  free  space. 

At  the  interface  both  the  tangential  component  of  electric  field  and  that 
of  magnetic  field  must  be  continuous.  For  the  magnetic  field,  Eqs.  (2.  1)  and 
(2.  4)  at  z  =  1  can  be  equated  to  obtain; 


£ 


a(k  )  e 

X 


-Jk  1 


"  e-j'^x^  dk 


x>  x) ^'1 


-Jk  X  , 
e  dk 


(2.  9) 


Similarly,  for  the  tangential  component  of  electric  field  Eqs.  (2.  2)  and  (2.  5) 
can  be  equated  and  thus. 


k  -ik  1  -ik  X 

a(k  )  e  e  ^  dk 

0 


(2.  10) 


c(k  e 


-jk 


£  -ik  X  ,,  ' 
e  ■’  X  dk 

X 


It  can  be  shown  that,  in  order  for  Eqs.  (2.  9)  and  (2.  10)  to  hold,  the 

I 

quantity  k  must  be  equal  to  k  .  This  result  is  not  surprising  because  it 
shows  that  the  phase  must  be  continuous  at  the  boundary,  which  is  also  the 
requirement  from  Snell's  law.  Since  the  integral  now  has  the  same  variable, 
the  quantities  inside  the  integral  can  be  equated,  noting  that: 
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k  =  k 

X  X 


^k  2  .  k  2 
P  X 


=  'Jk 


-  k 


0  X 

Thus,  from  Eqs.  (2.  9)  and  (2.  10)  results: 

-jk'  1 


-jk  1 
a(k^)  e  ^ 


r  jk  1 

b(k  e  ^  +  c(k  e  ^ 


(2.  11) 


and 


a(k  )  e 


-jk  1 


0 


P  L 


b<k  ^)  e 


-jk  1  ,  jk  1 

-  c(k  ^)  e 


]■ 


(2.  12) 


Equations  (2.8),  (2.  11),  and  (2.  12)  give  three  equations  relating  the 


three  unknowns.  In  solving  the  equations  for  st(k^), 


a(k^)  = 


27rE„  a  sin  k 

0  X  2 


jk  1 

..  z 


^z"0 

cos  k  1  +  j~k— J —  sin  k  ^1 
^  z  p  ^ 


(2.  13) 


Since  k  ,  k  ,  and  k  can  be  expressed  in  terms  of  k  .  the  right-hand  side 

X  z  Z  X 

of  Eq.  (2.  13)  can  be  expressed  completely  in  terms  of  k  • 

The  exact  solution  in  integral  form  has  now  been  obtained  [  see  Eq. 

(2.  1)]  .  Evaluation  of  this  integral  for  large  values  of  x  and  z  by  the  saddle- 
point  method  is  used  to  obtain  the  asymptotic  expression  for  the  far-field 
radiation  pattern. 

2.  2  Saddle-Point  Method 

The  saddle-point  integration  method,  also  known  as  the  method  of 
steepest  descent,  is  a  technique  to  obtain  an  asymptotic  approximation  for 
large  values  of  z  to  complex  integrals  of  the  type 

r  g(s)e^f<">ds 


where  g(s)  and  f(s)  are  independent  of  z. 

In  this  problem,  it  is  more  convenient  to  change  from  the  rectangular 
coordinate  system  to  cylindrical  coordinates.  Thus, 

k  =  k  sin  w 
X  0  ^ 


k  =  k„  cos  tp 
z  0  ^ 


dk 

dk  = 

X  dtp 


dtp  =  kQ  cos  ip  dtp 
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X  =  r  sin  6 

X  =  r  cos  0 

a(k  )  =  a(k.  sin  <fi) 
X  0 


where  (p ,  the  new  variable  of  integrations,  refers  to  the  direction  of  wave 
propagation,  whereas  0  refers  to  the  angle  of  observation  point.  Although 
both  0  and  (p  are  measured  from  the  z  axis,  0  is  independent  of  ip  just  as 
X  is  independent  of  k^  in  the  integration. 

Substitute  these  quantities  into  the  integral 

p  “  -j(k  X  +  k  z) 

H  =  \  a(k  )  e  ^  ^  dk 

y  J-=o  X  X 

and  obtain 


H 

y 


a(kQ  sin  (p)  k^  cos  tp 


-jk  r(cos  <i>  cos  0  +  sin  (p  sin  0) 
e  ’  dip 


=  ^  a(k  sin  (p)  k  cos  ip  e  ^  _  (2.  14) 

As  the  original  variable  of  integration  k^  goes  from  minus  infinity  to  plus 
infinity;  the  new  variable  <p  takes  the  path  C  shown  in  Figure  4.  The  shaded 
area  in  Figure  4  is  the  mapping  of  the  proper  sheet  of  the  k^  plane  into  the 
complex  <p  plane.  When  mapping  the  k^^  plane  into  the  (p  plane,  it  must  be- 
noted  that  the  imaginary  part  of 

k  =  k.  cos  0 
z  0 

must  be  negative  for  the  reason  of  convergence.  This  condition  is  satisfied 
in  the  shaded  region. 

The  next  step  of  this  procedure  is  to  distort  the  path  of  integration.  If 
the  path  C  is  continuously  deformed  into  another  path  B  in  such  a  way  that 
it  does  not  sweep  over  a  singularity  of  the  integrand,  then 

Ic  ■  Ib  . 

In  order  for  the  integral  to  converge,  it  is  also  required  that  the  two  ends  of 
the  path  stay  within  the  shaded  region. 

The  object  of  the  saddle-point  method  is  to  choose  a  path  of  integration 
so  that  the  quantity 

exp  |^R^(-jkQr  cos{0-<p} 

has  a  maximum  or  more  precisely  a  saddle-point,  and  falls  as  rapidly  as 
possible  on  either  side  of  this  maximum.  Also  along  the  path  near  the  maxi- 
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mum  the  imaginary  part  of  "j^Qr  cos(6-<f>)  is  to  be  constant,  or  in  other  words, 
the  phase  of  the  quantity 

g-jk^r  cos(e-(p) 

is  to  be  constant  along  the  path.  This  second  condition  insures  that  the  inte¬ 
grand  does  not  fluctuate  rapidly  between  positive  and  negative  values. 


Figure  4.  Complex  <p  Plane 
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The  saddle -point  occurs  when 


-^[-jV  cos(9-<p)]  =0 

or,  in  other  words,  at  =  0. 

If 

=  <^’l  +  jV’2  ' 

then 

-jk^r  cosIG-^^-jcp^)  =  k^r  ^sinh  sin(0-(;Pj^)  -j  cosh  cos(0-(^!j^)j  .(2.  15) 

The  imaginary  part  of  Eq.  (2.  15)  is  to  be  constant  along  the  new  path,  and  the  con¬ 
stant  can  be  evaluated  at  the  saddle-point  (p ^  =  0,  ^2~  gives  the  equation 

for  the  path 

cosh  <^2  cos(0-<p  j^)  =  1  .  (2.  16) 


It  should  be  noted  that  in  the  range  -T<<(p^<Tt,  there  is  more  than  one  possible 
path  which  satisfies  Eq.  (2.  16),  but  only  one  stays  within  the  allowed  region  in  the 
complex  (p  plane.  From  Eq.  (2.  16),  therefore. 


and 


cos(0-<;Pj^)  =  (cosh  <p^  ^ 

sin(0-</Jj^)  =  ±  fl  -  cos^(0-<;)  j^)j *  =  ±  tanh  (p^- 


(2.  17a) 


(2.  17b) 


The  minus  sign  must  be  taken  in  front  of  tanh  <p^  in  order  to  stay  on  the  allowed 
region. 

If  Eqs.  (2.  16)  and  (2.  17b)  are  substituted  into  Eq.  (2.  15),  then. 


-jk^r  cos(0-</j)  =  -k^r 


sinh^(/) 
cosh  tp^ 


(2.  18) 


Now  there  is  a  new  variable  (p^  which  goes  from  minus  infinity  to  infinity. 
It  can  be  seen  that  if  r  is  very  large,  the  quantity 


exp 


sinh  <p^ 
cosh</)2 


has  a  maximum  at  <^2  "  ^nd  decreases  very  rapidly  on  either  side  of  the  saddle- 
point.  Thus,  the  major  contribution  to  the  original  integral  comes  from  a  small 
region  of  the  path  near  the  saddle-point.  Therefore,  the  following  approximations 
near  the  saddle-point  may  be  made: 
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sinh^ij)^ 
cosh  (^2  "  ^2 

and 


d(jj 


dtp 

dtp 


dV’j  + 


3(^2 


d(p 
dtp^ 

-[i-i] 


dtp 


8(/j 

^<Po 


d(p , 


d(p 

3<P, 


d(/). 


dtp^  =  4J  dtp^ 


By  inserting  these  relations  into  Eqs.  (2.  18)  and  (2.  14),  it  can  be  seen  that 
the  major  portion  of  the  integral  is  given  by 


H 


j5r/4  ,  .  . 

^  m2  e  \  sin  tp)  k  cos  tp  e 

-n  ^  ^ 


■jV  •V^2' 


dtp^  (2.  19) 


where  u  is  a  point  on  the  path  near  enough  to  the  saddle-point  so  that 
2 


sinh^u 


cosh  u 


is  a  valid  assumption. 

If  r  is  made  very  large,  the  region  -u<tp^<M  can  be  made  very  small 
and  still  be  a  good  approximation  to  the  integral.  If  the  region  is 

small,  the  following  further  approximation  within  the  region  can  be  made: 

aikpSin  tp)  cos  tp  =*  aik^sin  9)cos  9 
The  quantity 

aikp  sin  9)  cos  9 

can  now  be  taken  outside  the  integral  and  the  integral  reduces  to 

_  iir/4  -k  rto  ^ 

Hy=>'v2  e-*  '  kg  cos  9  aik^sin  9)  e  j  e  0  ^2  dtp^.  (2.20) 


It  can  be  seen  that  if  the  limits  of  integration  are  extended  from  minus 
infinity  to  plus  infinity,  the  value  of  the  integral  changes  very  little  because 
r  is  assumed  to  be  very  large.  Thus  the  integral  is  approximated  further 
to  the  following: 

H  =  -/r  e^ 

y 


-jk  r  -k  (?  ^r 

Q  cos  9  aikp  sin  9)  e  J  e  ‘  ^‘(’2' 


The  integral  above  is  a  well-known  definite  integral  and  thus  the  approximate 
solution  is 
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H 


k„r  ^ k„  cos  0 


Q.  ^  cos  «  a(kQ  sin  0)  e 


“ik„r 
^  0 


where 


a(k|^  sin  0)  = 


27rEQ  a  sin(— 2 —  sin  0) 

V 

kp  cos  0  “2 —  sin  0 


jkyl cos  0 
e _ 

Enk' 

.  0  z _ 


(2 


cos  k  1  +  i  ,  „ 

z  ■’e  k_  cos  0 

L  p  0 


sin  k  1 
z 


.  22) 


I 


k  ^-k^^  sin^0 


In  the  analj'sis,  the  effect  of  poles  of  integrand  in  the  complex  (p  plane  has 
been  disregarded.  The  contribution  from  the  poles  can  be  included  in  the  usual 
marjier,  but  can  be  neglected  in  the  far-field  analysis  of  r  approaching  infinity 
unless  the  pole  lies  on  the  real  axis  of  the  <p  plane.  In  this  problem  there  was 
no  pole  on  the  real  axis,  and  thus  the  analysis  is  valid. 

The  expression  for  the  electric  field  can  be  obtained  in  a  similar  manner. 

If  and  E^  are  combined  to  obtain  a  single  component  Eg,  then 


E 


0 


0 

^  H 
0  y 


(2.  23) 


2.  3  Power  Pattern 


The  power  pattern  can  now  be  obtained  from  the  relation 


For  the  case  of  lossy  plasma  layer,  the  expression  for  the  power  pattern  be¬ 
comes  quite  complicated.  A  useful  form  of  power  pattern  cannot*  be  obtained  for  such 
a  case,  and  thus  only  the  expression  for  lossless  plasma  will  be  given  here.  However, 
the  equations  for  the  lossy  plasma  cases  are  presented  in  the  Appendix. 

I 

Since  the  quantity  k  ^  in  Eq.  (2.  22)  takes  on  real  and  imaginary  values 
for  the  lossless  case,  it  is  convenient  to  give  two  expressions  for  power 

I  I 

pattern;  one  valid  when  k  ^  is  real,  and  the  other  valid  when  k  ^  is  imaginary. 

Let 

1  '  ,  2  .  2  _ 

k  =  vk  -  k„  sin  0 
zr  p  0 

Y 


where 
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Then 


and 


k  =  k  for  cos  6  >_£- 
z  zr  u) 


=  ik  .  for  cos  6  <  — £ 


It  should  be  noted  that  when  the  plasma  layer  is  lossy,  the  quantity 
k  is  complex  instead  of  pure  real  or  pure  imaginary. 

^  I 

By  substituting  the  relations  above  for  k  ^  into  Eq.  (2.  22),  the  following 
expressions  for  the  normalized  power  pattern  are  obtained; 


P  =  — 
n  r 


II  sitff 

‘J 


sin  6 


2  ' 

cos  k  1 + 
zr 


L 


2  '  2 
e.'^k 
0  zr 

2,  2  2  . 
£  k„  cos  e 
P  0 


•  2  , 
sin  k 


zr 


(2.  24a) 


for 


and 


for 


cos  0  >— 

(jC 


p  _  1^,  sinV  2 


Sin  e 


koa 


sin  0 


2  r 


cosh  k  .1  + 
zi 


e  2  k'  2 

0  ^  zi 

e  ^  k„^  cos^e 

P  0 


sinh  k  21^ 


(2.  24b) 


cos 


0  <-£ 


In  normalizing  the  power  pattern,  it  was  assumed  that  the  voltage  across  the 

gap  is  constant  for  all  cases. 

2 

Since  cos  0  is  never  greater  than  unity,  only  Eq.  (2.  24b)  is  needed 
when  ^p  is  greater  than  unity.  However,  when  ‘^p  1  there  is  a  definite 

angle  at  which  k  ^  changes  from  real  to  imaginary  quantity.  This  angle  will 
be  called  the  cutoff  angle.  The  cutoff  angle  is  given  by 


0  =  cos 

c 


1  ^ 
■1_2_ 


In  Figures  5  and  6, some  radiation  patterns  are  plotted.  In  plotting  them,  it 
was  assumed  that  the  slit  was  infinitesimally  narrow.  This  assumption 
gives  the  relation 

— sin  9 
sin  2  , 


sin  0 


I 
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O) 

In  Figure  5,  pKSwer  patterns  are  shown  for  various  values  of  with 

fixed  plasma  thickness.  The  cutoff  angles  are  indicated  on  the  graph.  An 

interesting  feature  Oi  the  patterns  is  that  when  ^p  <  1,  the  patterns  have 

oi 

peaks  just  inside  the  cutoff  angle.  For  the  values  of  9  greater  than  the  cut¬ 
off  angle,  the  amplitudes  decrease  rapidly  with  increase  in  the  angle. 

VVhen-^  ^  1,  the  amplitudes  of  the  pattern  are  greatly  attenuated.  This 
effect  is  due  to  the  propagation  constant  in  the  plasma  given  by  the  expres¬ 
sion 


which  becomes  imaginary.  Thus  the  waves  are  attenuated  exponentially 
within  the  plasma  layer. 
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Figure  6  shows  the  power  patterns  plotted  as  functions  of  plasma  layer 
thickness.  For  the  particular  value  of  chosen,  the  cutoff  angle  was  37°. 

O) 

Since  the  power  patterns  are  symmetric  in  9,  only  half  of  each  pattern  is 
plotted.  The  peaks  fall  approximately  at  the  same  angle  and  they  are 
higher  for  thicker  plasma  layers. 


Figure  6.  Power  Patterns  for  Values  of  k^l  (TM  Wave  Case) 


Since  the  ground  plane  is  never  infinite  in  practice,  an  important 

cu 

relation  is  the  attenuation  as  a  function  of _ P  for  9  =  0.  The  plot  of  attenua- 

U)  “ 

tion  versus _ E  with  the  plasma  thickness  as  a  parameter  is  given  in  Figure 

u) 

7. 


20 


ATTENUATION  (DB.) 


2.4  Plane  Wave  Model  and  Transmission  Line  Analogy 

To  get  some  physical  feeling  for  the  problem,  consider  the  following 
system. 

Let  a  fictitious  plane  wave  source  at  z  -  0  which  constrains  the  x  com- 
pwnent  of  the  electric  field  to  be  of  the  form 

-jk_x  sin  9 

e^  (z  =  0)  =  Eq  cos  9  e  , 

where  =  wn/pqCq  .  (Small  letters  are  used  to  denote  the  field  compo¬ 
nents  for  the  plane  wave  model.  ) 
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Figure  8.  Plane  Wave  Source 


If  freedom  to  choose  0  is  assumed,  and  the  source  extends  from 
X  =  -oo  to  X  =  “,  then  such  a  source  propagates  a  plane  wave  in  the  6 
direction  with  the  magnetic  field  given  byt 

ro  -jk^  (x  sin  6  +  z  cos  9) 

H  =  E  e  I  .  (2.  25) 

0  u  y 

Thus  for  any  6  the  choice  is  between  -  w/2  and  ff/2;  the  source  propagates 
a  plane  wave  with  the  electric  field  amplitudes  equal  to  E^. 

Now  suppose  that  a  slab  of  plasma  is  placed  against  the  source  as 
shown  in  Figure  9. 


iZ 

II 

77777777777777777777777, 

777777777m 
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PLANE  WAVE  SOURCE 


Figure  9.  Plane  Wave  Source  with  a  Plasma  Slab 

To  solve  for  the  fields,  assume  plane  wave  solutions  in  the  two  regions 
of  interest.  Inside  the  plasma  layer  there  is  a  solution  of  the  following  form; 
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t  1 

h  =  H 


y  0  LT 


-jk  +  r  e  J  e 


1  9h  H  „k 

JL  =  0  z 


X 


-ik  z  ik  z 
z  T-  z 
e  -re 


-ik  X 
X 


'  _  1  ^  ^  -h'o^'x  I  z"  ^ 

®  z  jcje  9x  oie  ®  ^ 

•'  P  P 


-ik  X 

•'  X 


where 


k  =  k  sin  0 

X  p 


k  =  k  cos  0 
z  p 


In  the  region  z  >  1. 


h  =  e 

y  0 


■j(k  z  !  k  x) 
z  X 


H„k  -i(k  z  +  k  x) 
e  =— S-2-  e  ^ 

X  we. 


■^0*^x  -.i(k  z  +  k  x) 

e  - -  z  X 

z  wCq  e 


where 


k  =  k„  sin  0 
X  0 


k  =  k„  cos  0 
z  0 


At  z  =  0.  the  following  relation  from  the  txDundary  condition  on  the  x 
compxjnent  of  the  electric  field  results; 


I  I  1 

-ik„x  sin  0  H  .k  -ik  x 

TT.  q  O  0z,,„,  ■’  X 

E„  cos  0  e  =  -  [  1  -  r]  e 


0 


W€ 


Since  the  relation  must  hold  for  any  x.  then 


k.  sin  0  =  k  =  k  sin  0 

0  X  p 


and 


I  ? 

H  „k 
0  z 


[  1  -  r)  . 


(2,  26a) 

(2.  26b) 

(2,  26c) 


(2,  27a) 

(2.  27b) 

(2  27c) 


Eq  cos  0  = 


(2.  28) 
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From  the  boundary  conditions  at  z  =  1  ,  the  following  relations  are 
obtained; 


and 


-jk  1  , 

Hq  "  =  «  0 


-jk  1  jk  1 
e  ^  +  r  e 


k  -jk  1  h'  k' 

T,  z  z  _  0  z 

—  e  - 


0  £ue 


0 


we 


P  L 


c  -re 


(2.  29) 


(2.  30) 


There  are  three  equations  (2.  28),  (2.  29),  and  (2.  30),  and  three  unknowns 
H  Q,  Hq  and  T.  Solving  for 


H„  .  E„  JS.  ibl 
0  0  fi.  e 


■  z^O 

cos  k  1  +  j-r— r —  sin  k  1 
z  k  c  z 

z  p 


(2.  31a) 


Thus  the  amplitude  of  the  magnetic  field  in  the  free  space  region  is  given  by: 


0  ^0 

: - - = 

2  ,  '  ,  ,  z  ^0  .  2  ,  '  , 

COS  k  1+ - K — TT  sin  k  1 

z  k  2  z 

z  p 


(2.  31b) 


The  interesting  part  of  this  result  is  that  as  far  as  the  9  dependence  is 
concerned,  the  expression  for  is  exactly  the  same  as  the  far-field 
pattern  obtained  for  the  narrow  slit  case  [  see  Eq.  (2.  22)]  .  The  similarity 
of  the  results  suggests  interpretations  for  some  features  of  the  p)ower  pat¬ 
terns  in  terms  of  the  plane  wave  model.  However,  in  order  to  make  the 
interpretations  easier,  analogies  can  be  drawn  with  the  transmission  line. 

The  transmission  line  analogue  of  the  plane  wave  model  is  found  to  be 
as  shown  in  Figure  10.  (See  Reference  4  for  transmission  line  analogues 
of  plane  waves.  ) 


Figure  10.  Transmission  Line  Analogue 
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The  analogue  gives  the  z  variation  of  the  field  components  e  ,  h  , 

X  y 

e^  and  h^.  If  the  system  is  solved  for  h^,  then 


rrr- 


h  =  E  - 

y  °  >^0  r- 


-ik^(z-l ) 


,  '  ,  ,  z^O  .  ,  '  , 

cos  k  1  +  — - —  sin  k  1 
z  •’  k  e  z 

z  P 


The  power  propagated  per  unit  area  can  be  computed  from 


P  =  2 


When  e  =  0,  then 


J3. 


Let  us  investigate  the  peaking  effect  observed  in  Figures  6  and  7.  It 
should  be  noted  that  these  peaks  are  greater  in  amplitude  than  the  case  with¬ 
out  the  plasma  slab. 

From  Eq.  (2.  31), 


h  =  E 


Op 


1 


1  -  1  - 


<2  2 
^z  ^0 

k  2 

z  p 


Thus  |h  I  can  become  greater  than  E^^^— ^ 

y  ®  0  Pq 

when 


2  ' 

Isin  k  ^1 


k  e 

0  <  ^  ■  <  1  . 
k  € 
z  p 

This  condition  is  indeed  satisfied  in  the  range 


.  -1 
sin 


e  <  e. 

p  0 


and  within  this  range  of  angles,  the  power  radiated  is  greater  than  the  "no 
plasma"  case. 

Similar  relations  are  obtained  for  H  in  the  narrow  slit  problem  [  see 
Eq.  (2.  22)]  . 

In  terms  of  the  transmission  line  analogue,  the  condition 


k  €„ 
z  0 

k  e 
z  p 


<  1 
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is  analogous  to 


or 


I 


By  adding  a  transmission  line  of  lower  characteristic  impedance  to  a  matched 
line,  the  impedance  seen  from  the  source  is  always  reduced.  Hence,  more 
power  is  delivered  to  the  load. 

The  transmission  line  analogy  also  shows  why  fluctuations  in  the  power 

pattern  for  thick  plasma  layers  are  observed.  As  0  changes  from  0  to  0  , 

I  c 

the  propagation  constant  k  ^  in  the  transmission  line  analogue  changes  from 

kp  to  0.  Thus  when  k^l  »  1,  the  length  1  changes  from  many  wavelengths  to 
zero  wavelength  as  0  is  changed  from  0  to  0  .  In  other  words,  the  impedance 
seen  from  the  source  fluctuates  as  0  is  varied  and  fluctuations  in  the  power 
pattern  result. 

CO 

When  the  angle  0  is  greater  than  0^  or  when  ^  1,  the  transmission 
line  analogy  no  longer  holds.  In  such  a  case,  the  plane  waves  within  the 
plasma  layer  are  non-uniform  and  are  attenuated  in  the  z  direction. 


3.  TRANSVERSE  ELECTRIC  WAVE 


3. 1  Solution 

Now  consider  the  case  of  the  electric  field  polarized  along  the  gap  in  y 
direction.  The  geometry  for  this  problem  is  shown  in  Figure  11.  The  field- 
strength  distribution  across  the  gap  is  assumed  to  be  co-sinusoidal  as  shown 
in  Figure  12.  This  particular  configuration  corresponds  to  a  waveguide 
opening  into  half  space  with  the  waveguide  excited  by  TE^^^  mode. 


Figure  11.  Radiative  Geometry  for  the  TE  Wave  Case 


P'igure  12.  Electric  Field  Amplitude  Distribution  at  z  =  0  for  the 
TE  Wave  Case 
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The  method  used  to  solve  this  problem  is  very  similar  to  the  one  used 
in  the  previous  problem.  In  this  particular  polarization  and  geometry,  it 
can  be  seen  that 


1  1  t 

E  =  E  =  H  =  0 
X  z  y 


There  is  only  one  component  of  electric  field  and  in  the  free  space  it  is 
represented  by  the  following  integral; 

-j(kx+k  z) 

E  (x,  z)  =  \  a(k  )  e  ^  ^  dk  (3.  1) 

y  J-oo  X  ^ 

The  solution  above  is  identical  to  the  one  assumed  in  the  previous  problem 
except  this  time  the  solution  is  assumed  for  the  electric  instead  of  magnetic 
field.  The  corresponding  magnetic  field  comfwnents  in  free  space  region  are 
given  by  the  following  relations: 


H  (x 

X 


^  3E  p*  k 


-j(k  X  +  k  z) 

a(k  )  e 


dk 


(3.2) 


and 


.  8E  po 

“jwpQ  Sx  '  J_ 


>00  k  -3(k^x  +  k^z) 

a(k  )  e 


*  dkj^  .  (3.  3) 


Inside  the  plasma  medium,  the  solution  is  represented  in  terms  of  out¬ 
going  waves  and  reflected  waves; 


E  y  (x,  z) 


r. 


-jk  ,z  ,  jk  1  -jk  X  , 

e  dk  ^  (3.4) 


^)  e  ^  +  c(k  ^)  e  ^ 


and  the  corresponding  magnetic  field  is  given  by; 

8e' 


H  (x,  z)  =-:^ -  — 5-- 

X  jwPq  9z 


-^^[b(k' 


.)  e 


-jk  z 


-  c(k 


X  J 


-jk  X  , 
e  dk 


X 

(3.5) 
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H  (x,  z) 
z 


-1 


r  ^ 

J.OC 


+  c(k 


jk 


(3.  6) 


At  the  ground  plane  the  tangential  component  of  electric  field,  E  (x,  0), 
must  be  0  everywhere  except  at  the  gap  where  the  amplitude  distribution 
becomes  co-sinusoidal  as  discussed  earlier.  Equation  (3.  4)  for  the  electric 
field  becomes  at  the  ground  plane; 


E  y(x,  0) 


b(k  )  + 


-jk  X  , 

e  ^  dk 


which  indicates  that  the  quantity 


(3.  7) 


is  the  Fourier  transform  of  the  co-sinusoidal  electric  field  amplitude  dis¬ 
tribution.  By  equating  the  Fourier  transform  to  the  expression  above,  then 


|b<k  +  c(k 


]  = 


2^Eq  cosk^^ 

{if  - 


(3.  8) 


Equation  (3.  8)  gives  one  of  the  three  required  relations  to  solve  for  a(k  ). 

At  the  interface  between  the  plasma  and  the  free  space  continuity, 
conditions  on  the  tangential  components  of  electric  and  magnetic  fields 
result.  Thus,  the  two  electric  fields  E  1)  and  E^(x,  1)  are  equated  and 
also  the  two  magnetic  fields  H  (x,  1)  and  H  (x,  1).  In  equating  them,  it 
can  be  seen  again  that  k  ^  must  be  equal  to  k^.  The  equality  provides  the 
equation  of  the  quantities  inside  the  integral  as  in  the  previous  problem. 

By  equating  Eqs.  (3.  1)  and  (3.  4)  at  z  =  1,  then 


-jk  1  ,  -jk  1  1  jk  1 

a(k  )  e  ^  =  b(k  )  e  ^  +  c(k  )  e 

X  X  X 


(3.  9) 


for  the  electric  field  continuity.  For  the  magnetic  field  continuity  relation, 
from  Eqs.  (3.  2)  and  (3.  5)  results 


a(kx) 


=  k  lb(k  )  e 


-jk  1 
•’  z 


-  c(k  „)  e 


jk 


1 


(3.  10) 


Equations  (3.  8),  (3.  9)  and  (3.  10)  give  three  equations  and  three  unknowns 
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enabling  us  to  solve  for  a-(k^).  By  solving  the  equations  for  there 

results 


27rE„  cos  k 


X  2 


cos  k  1  +  j- 
z  % 


sin  k  1 
z 


(3.  11) 


I  I 

As  before,  Snell's  law  allows  the  quantities  k  k^  and  k  ^  to  be  expressed 
in  terms  of  k  . 

X 

The  far--field  approximate  expression  for  the  electric  field  can  be  ob- 
tcdned  by  the  saddle -point  method,  and  the  result  is  as  follows: 


E  “  k  cos  0  a(k  sin  0)  e  ^  (3.  12) 

y  K^r  u  U 

where 

kQr  »  1. 

If  Eq.  (3.  11)  is  used  to  evaluate  alk^^  sin  0),  the  following  substitutions  can 
be  used: 

I 

k  =  k  =  k.  sin  0 
X  X  0 

k  =  k.  cos  0 
z  0 

»]  *  . 


3.2  Results 

The  normalized  power  pattern  can  be  obtained  from  the  relation 


Again,  only  the  lossless  case  is  considered.  In  the  lossless  case,  the 
1 

quantity  k  ^  takes  on  either  real  or  imaginary  values  as  in  the  case  of  the 
transverse  magnetic  wave  studied  in  the  previous  section.  Thus  two 
expressions  are  given  for  the  power  pattern. 

Let 

i 

k  =  fk  ^  -  k_^  sin^  ol  for  k  ^  k.  sin  0 

zr  [_  p  0  J  P  0 
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and 

k  .  =  k 
zi  L 

0  ® 

-  k 

1 

1 

where 

II 

r  ^  "1 

^  a 

P  0  ! 

2  1 

21  ^ 


for  k  ^  sin  G 
P  0 


Then  the  following  expressions  for  the  power  pattern  can  be  obtained  from 
Eqs.  (3.  11)  and  (3.  12).  The  power  is  normalized  so  that  the  electric  field 
amplitude  between  the  slit  is  kept  constant  for  all  cases.  Furthermore, 
the  value  of  i  is  given  to  the  power  radiated  in  0  =  0  direction  for  the  "no 
plasma"  case. 


for 


n  r  \  a/ 


k  >  k„  sin  G 
P  0 


for 


cos\  2  sin  G  / 

2 

r 

2„ 
cos  G 

/ir\  ^  1  2  .2 

Lu)  •‘^0 

2  2 

o  1  k^  cos  0 

cos  k  1  +  — “  '1 — 

zr  '  ^ 

k 

^  zr 

2  ' 

sin  k  1 
zr 

(3.  13a) 

and 


n  r  \a7 


cos 


2  „ 
cos  9 


cosh  k 


^  zi  ' 


sinh  k  .1 
zi 

(3.  Ub)' 


k  <  k.  sin  0 
P  0 


If  the  plasma  layer  is  lossy,  it  is  necessary  to  use 

2 
u 

1  + 


k  =  k 
P  0 


u 

_ P  . 

(i'+jw)iw 


to  obtain  k 


In  such  a  case,  k 


^  - ,  ..  ^  is  complex  and  a  simple  expression  for 

power  cannot  be  obtained.  However,  the  exact  expressions. for  the  lossy  plasma 
case  are  given  in  the  Appendix. 

For  the  particular  electric  field  |X>larization  considered  in  this  section, 
radiation  from  an  infinitesimally  narrow  slot  is  not  possible.  Thus,  for 
the  calculation  of  the  radiation  pattern,  a  value  of  .  8  was  assigned  to  ^  , 
the  ratio  of  slot  width  to  free  space  wavelength.  The  value  assigned  is  a 
typical  value  for  the  ratio  of  waveguide  width  to  free  space  wavelength  for 
a  TEj^q  waveguide  mode. 
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The  radiation  patterns  are  shown  in  Figures  13  and  14.  Unlike  the 
transverse  magnetic  case  treated  in  the  previous  section,  the  power  pat¬ 
tern  without  the  plasma  layer  is  not  omni-directional. 


Figure  14.  Power  Patterns  for  Various  Values  of  k^l  (TE  Wave  Case) 
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Again,  analogies  with  the  transmission  line  can  be  drawn  as  was  done 
in  Section  2.  However,  the  characteristic  impedances  of  the  equivalent 
transmission  lines  must  be  re-defined.  For  instance,  it  is  necessary  to 

identify  with^^ —  instead  of  —  .  (See  Figure  10  and  the  subsequent 

‘^z  0 

analysis.  ) 

I 

Since  the  quantity  k^/k  ^  in  the  denominator  of  Eq.  (3.  13a)  is  always 

greater  than  unity  for  real  values  of  k  the  radiated  power  never  exceeds 

that  of  the  "no  plasma"  case.  As  a  result,  peaks  are  observed  whenever 
1 

k  1  =  nir  where  n  is  a  non-zero  integer.  Whenever  k  1  <  tt,  the  only  peak 
z  p 

is  the  one  at  0  =  0  and  the  power  pattern  decreases  monotonically  with  0. 
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4.  ANISOTROPIC  PLASMA 
4. 1  Introduction 

In  the  presence  of  a  magnetostatic  field,  the  plasma  medium  becomes 
8  9 

anisotropic.  '  In  general,  with  the  magnetostatic  field  in  the  z  direction 
in  a  right-hand  rectangular  coordinate  system,  the  relative  dielectric 
permittivity  can  be  described  by  the  following  tensor  (see  Section  1.  2). 


Because  the  off  diagonal  elements  of  the  tensor  are  complex  quantities, 
the  electromagnetic  waves  in  such  a  medium  propagate  in  a  very  complex 
mainner.  Thus  only  the  special  case  of  a  very  high  magnetostatic  field  is 
studied. 

When  the  magnetostatic  field  becomes  very  high,  the  dielectric  tensor 
approaches  the  following  expression; 


The  off  diagonal  elements  go  to  0  and  the  only  remaining  element  different 
from  the  free  space  permittivity  is  the  one  associated  with  the  direction  of 
the  magnetostatic  field  vector. 

In  the  analysis  to  follow,  it  is  assumed  that  the  magnetostatic  field  is 
strong  enough  so  that  the  above  tensor  is  a  valid  representation. 

4.  2  Propagation  in  Anisotropic  Medium 

At  this  point,  an  investigation  of  electromagnetic  wave  propagation  in 
an  anisotropic  medium  is  appropriate. 

Consider  the  geometry  shown  in  Figure  15.  The  space  is  filled  with 
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an  anisotropic  dielectric  which  can  be  described  by  a  dielectric  tensor  of 
the  following  form: 


0 

0 


0 

£ 

y 

0 


Assume  that  a  uniform  plane  wave  travels  in  the  0  direction  and  per¬ 
pendicular  to  the  y  axis,  and  that  the  magnetic  field  is  polarized  in  the  y 
direction. 

Since  the  medium  is  non-gyrotropic,  this  configuration  is  as  follows: 


H 


X 


=  0. 


From  the  equation 


V  X  H  , 

we  obtain  the  relations 
1  9H 

E  = 

ja:£^ 


and 


3z 


8H 


E  = 
z  3a)€ 


-IL 


ax 


The  other  curl  equation 
V  X  E  = 

yields  the  following  relation: 
-1 


H  = 

y 


aE 


az 


8E 


ax 


(4.  1) 


(4.  2) 


(4.  3) 


By  substituting  Eqs.  (4.  1)  and  (4.  2)  into  (4.  3),  the  following  partial  dif¬ 
ferential  equation  for  is  derived: 


H 

y 


(4.  4) 


The  general  solution  to  such  an  equation  is  given  by: 
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Equation  (4.  5)  is  the  relation  that  must  hold  between  k  and  k  in  the  mode 

X  z 

of  propagation  considered. 

The  propagation  constant  k,  as  a  function  of  6,  can  be  obtained  by 
substituting  k  =  k  sin  0  and  k  =  k  cos  0  into  Eq.  (4.  5).  It  is  found  that 

X  z 


e  €  w  u 

X  z  ^0 _ 

■  2  2 
I  e  cos  0  +  €  sin  0 
Lx  z 


Figure  15.  Plane  Wave  in  Anisotropic  Medium 
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4.  3  Solution 

The  problem  of  finding  the  radiation  pattern  with  the  anisotropic  plasma 
is  approached  in  much  the  same  way  as  in  the  previous  cases.  It  is  seen  that 
there  are  only  two  interesting  configurations.  In  both  cases,  the  polariza¬ 
tion  of  the  electric  field  at  the  gap  is  in  the  x  direction  as  in  the  problem  of 
Section  2.  In  one  case,  the  magnetostatic  field  vector  is  in  the  z  direction, 
and  in  the  other,  it  is  in  the  x  direction.  Any  other  combination  of  magneto¬ 
static  field  direction  and  electric  field  polarization  reduces  either  to  the 
case  with  no  plasma  or  to  the  case  already  studied  due  to  the  fact  that  two 
of  the  diagonal  elements  of  the  tensor  have  the  value  of  Thus,  only  the 
two  cases  mentioned  are  considered. 

The  geometry  and  the  excitation  of  this  problem  are  the  same  as  de¬ 
scribed  in  Figures  2  and  3.  Since  the  plasma  medium  is  non-gyrotropic,  it 
is  expected  that. 


H  =  H  =  E  =  0 
X  z  y 


I  I  I 

H  =  H  =  E  =  0 
X  z  y 


assuming  of  course  that  the  magnetostatic  field  points  either  in  x  or  z 
direction.  As  before,  the  solution  for  magnetic  field  in  free  space  is 
assumed  to  be  of  the  following  form; 


p"  -j(k  X  +  k  z) 

H  =  \  a(k  )  e  ^  ^  dk 

y  J  „  X  x 


(4.  6) 


with  the  corresponding  electric  field  components 
E  =  -\  — ^  a(k  )  e  ^  ^  dk 

and 


(4.  7) 


k  -MK  : 

=  r  ^  a(k  )  e  ^ 


-j(k^x  +  k^z) 


dk 


(4.  8) 


Inside  the  plasma  medium,  a  solution  identical  to  the  one  assumed  in 
Eq.  (2.  4)  of  Section  2  is  assumed  for  the  magnetic  field  as  shown; 


h'  -f 
y 


..  ^  -jk'  z  ,  jk  z'’ 

b(k  e  +  c(k  e 


-jk  x  , 

e  ^  dk  ^  .  (4.  9) 


I 
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However,  this  time  it  must  be  noted  that  there  is  a  different  relation  between 

I  I 

k  and  k  ,  namely; 

X  z  •' 


,  2 


k 


X 


2 

cj  H 


,  2 


0  X 


(4.  10) 


The  eleetric  field  in  the  plasma  medium  is  given  by; 


'  -  f 

X  J  ^ 

X 

I 

t  r*  k 
=  \ 

Z  J _ cje„ 


and 


.  -jk  z  ,  jk  z 

b(k  e  -  e(k  e 


1  "jk  z  ,  . 

b(k  )  e  ^  +  e(k  )  ^  z^ 
X  ■  X  e 


■jk  X  , 

j  dk  ^  (4.  11) 


-jk  X  , 

e  ^  dk  ^  .  (4.  12) 


The  proeedure  heneeforth  for  finding  a(k^)  is  exaetly  as  that  in  Seetion 
2,  and  thus  a  detailed  analysis  is  omitted.  The  Fourier  transform  teehnique 
is  used  to  mateh  the  boundary  eonditions.  As  before,  the  eontinuity  eon- 
ditions  on  the  eleetrie  and  the  magnetie  fields  yield  Snell's  law  giving  the 
relation,  k  =  k  . 

X  X 

The  resulting  funetion,  atk  ),  is  identieal  to  Eq.  (2.  13)  exeept  that  e 

X  p 

is  replaeed  by 


2jrE  a  sin  (k  -z) 

a(k  )  = - 2 — ; - ^ 

^  k  (k  -) 

•^z'  X  2' 


jk  1 
-  z 


'  ^  z^'O  ' 

eos  k  1  +  j  yT-  sin  k  1 

Z  X 


(4.  13) 


From  the  relation 
(4.  9)  to  get 


k 


z 


k  ^  can  be  evaluated  in  terms  of  k^  using  Eq. 


O^x 


1  - 


~2 - 


(4.  14) 


The  values  of  and  depend  on  the  direction  of  the  magnetostatic  field 

vector.  The  relation  between  k  and  k  is  given  by: 

X  z 


k 

z 


1 

9  2  * 

f'o'o  •  ‘^x  > 


(4.  15) 


Thus,  a(k^)  can  be  expressed  entirely  in  terms  of  k^. 

By  applying  the  saddle-pwint  method,  the  far-field  pattern  is  obtained 
and  is  as  shown  below: 
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H 

y 


.n/21 


cos  6  a(kQ  sin  G)  e 


-j<V  -  f) 


(4.  16) 


The  quantity  a(kQ  sin  0)  can  be  evaluated  using  Eq.  (4,  13)  and  the  auxiliary 
relations  in  Eqs.  (4.  14)  and  (4.  15). 


4.4  Results 

In  the  following  analysis,  it  is  assumed  that  the  plasma  layer  is  lossless. 
The  power  pattern  for  the  lossy  case  can  be  obtained  by  substituting  the  proper 
values  of  e  and  e  into  Eqs.  (4,  13),  (4.  14),  (4.  15),  and  (4.  16). 

X  z 

(a)  Magnetostatic  field  vector  in  the  z  direction 
When  the  magnetostatic  field  is  z  directed,  then 


By  substituting  these  quantities  into  Eq.  (4.  14)  and  setting  k^  =  k^  sin  G, 
there  results 

z 


k  =  k 
z  0 


cos  0  - 


(4.  17) 


The  expression  above  shows  that  when  — £  is  greater  than  unity,  k  is 

a.’  2 

always  real.  However,  when  1,  k  ^  takes  on  either  real  or  imaginary 

values.  Thus  it  is  again  convenient  to  obtain  two  expressions  for  the  power 
pattern. 


Let 


=  k 


u 

e 


(jJ 

1-^ 


for  cos  0 


u 


and 


k  .  =  k. 
zi  0 


-  cos  0 


CJ 


u 


for  cos  0 


By  using  these  relations  with  Eqs.  (4.  13)  to  (4.  16),the  following  expressions 
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for  the  normalized  power  pattern  are  obtained: 


kQa  n 

sin  sin  0) 

2 

1 

k  a 

— 5—  sin  0 

L  2  j 

2  '  / 
cosh  k  .1  + 

L  \ 

1  2 
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kocos0j  k^.lj 
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amd 


P  =  - 
n  r 


cos  0  <  — <  1 
u 


sin  sin  0) 


kQa 


sin  0 


cos  k  1 + 
zr 


zr 


k^  cos  0 


sin  k  1 
zr 


(4.  19) 


for 


k  =  real  number, 
z 


Equations  (4.  13)  and  (4.  19)  give  the  power  pattern  for  the  case  of 
infinite  magnetostatic  field  directed  parallel  to  the  z  axis. 

The  radiation  patterns  as  function  of  are  shown  in  Figure  5.  It  is 
assumed  that  the  slot  is  infinitesimally  narrow.  Qualitatively  the  patterns 
are  quite  similar  to  those  in  Figure  5,  the  case  of  no  magnetic  field.  How¬ 
ever,  unlike  the  patterns  of  Figure  5,  the  peaks  are  higher  for  larger 
values  of  — i-  and  the  normalized  power  is  unity  in  the  direction  of  the 
magnetostatic  field  even  when—i  is  greater  than  unity.  This  latter  effect 
is  expected  because  the  propagation  constant  in  the  direction  of  magneto¬ 
static  field  reduces  to  that  of  the  free  space.  The  wave  traveling  in  this  direc¬ 
tion  has  its  electric  field  polarized  perpendicular  to  the  direction  of  the  mag¬ 
netostatic  field.  Since  the  dielectric  permittivity  associated  with  the  direction 
perpendicular  to  the  magnetostatic  field  is  the  wave  propagates  with  the 
free  space  propagation  constant. 

(b)  Magnetostatic  field  vector  parallel  to  the  x  axis 
Now  consider  the  case  of  magnetostatic  field  parallel  to  the  x  axis.  In 
this  case,  2 


€ 

Z 


€ 


0 


By  substituting  these  relations  into  Eq.  (4.  14)  and  setting  sin  9 
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1 
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k  =  k.  cos  0 
z  0 

L 

Clearly,  k  j.  real  when <  1  and  is  imaginary  when-^>  1  . 


Let 


k  =  k„  cos  0 
zr  0 


1  -- 


for  — ^  <  1 

O) 


and 


k  .  =  k„ 
zi  0 


W 

^  - 


cos  0 


for  —2  >  1  . 
oi 


U) 

Figure  16.  Power  Patterns  for  Various  Values  of — E.  with  the 
Magnetostatic  Field  in  the  z  Direction 


Then  from  Eqs.  (4.  13)  to  (4.  16)  the  following  expressions  for  the  normalized 
power  pattern  are  obtEiined: 
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cosh^  k  1  +  / - - - ^  \  sinh  k  .1 

zi  /  ^  2  \  zi 

u 


for 


1  . 


Shown  in  Figure  17  are  the  patterns  with  the  magnetostatic  field 
parallel  to  the  x  axis.  Unlike  the  previous  cases  studied,  power  radiation 
is  greatest  along  the  ground  plane.  In  addition,  the  amplitude  never  ex¬ 
ceeds  unity.  As  expected,  the  normalized  power  radiation  is  unity  in  the 
direction  of  the  magnetostatic  field. 


O) 

Figure  17.  Power  Patterns  for  Various  Values  of — Ewith  the 
Magnetostatic  Field  in  the  x  Direction 


(4.  20a) 


(4.  20b) 
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Appendix 


Radiation  Patterns  for  Plasma  Layers  with  Losses 
(No  Static  Magnetic  Fields) 

The  jx)wer  patterns  for  the  lossy  plasma  layer  can  be  derived  from  the 
general  expressions  for  the  transverse  magnetic  and  electric  wave  cases 
[  Eqs.  (2.  22)  and  (3.  12)  respectively)  by  allowing to  become  complex 
in  accordance  with  Eq.  (1.  8).  The  following  results  are  obtained  after  some 
algebraic  manipulation: 

(A)  Transverse  Magnetic  Wave 


(Al) 


P  =  - 
n  r 


sin  (kg  sin  0) 


ko  ^  sin  e 


4(A^  +  B^) 

C  +  D  +  E  +  F 


2 

where  =  n  =  A  -  jB 


C  =  e 


2al 


D  =  e 


-2al 


^(Akp  cos  6  +  +  (Bkp  cos  0  +  a)^^ — ^ 

(-Bk^  cos  e  +  a)^  +  (Akp  cos  0  -  ^ 


E  =  4  sin  201 


Q0  - 


AB  kp^  cos^  0 


>1- 

J  kJ 
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F  =  2  cos  2^1 


1, 

a  =  ^(A  -  sin^  6)^  +  sin  y 


^  cos^  e  +  cos^  e  -  -  /3^ 


[(A  -  sin^  e)^  +  B^  ]  ^ 


(3  =  k^  -  sin 


1  ,  -1  /  B 

7  =  2  tan  - 


cos  7 


2 

A  -  sin  e 


(B)  Transverse  Electric  Wave 


(A2) 


=  -  (^) 
n  \  a  / 


cos  (kp  I  sin  e) 


2  •  2  . 
sin  8 


4  cos^  6  [ 

[  C  +  D  +  E  +  F] 


where 


P  =  n^  =  A  -  jB 


2al 


[« 


C  =  e  I  (;3  +  kp  cos  8)  + 


D  =  e 


■2al 


^(;3  -  kp  cos  8)^  + 


E  =  2(^^  +  -  kp^  cos^  8)  cos  2j31 

F  =  -4  a  kp  cos  8  sin  2^1 


sin  7 


a  =  kp  [(A  -  sin^  8)^  +  B^] 

i 

^  =  kp  ^(A  -  sin^  8)^  +  B^J  cos  y 


1  .  - 1 
7  =  2  tan 


B 


A  -  sin  8 


It  should  be  noted  that,  since  B  is  a  positive  number  in  the  notations  used, 
the  following  relation  for  y  holds: 
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